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We investigate the groundstate properties of a recently proposed model for a topological Kondo 
Insulator In one dimension (i.e., the p-wave Kondo-Heisenberg lattice model) by means of the Den¬ 
sity Matrix Renormalization Group method. The non-standard Kondo interaction in this model 
Is different from the usual (i.e., local) Kondo interaction in that the localized spins couple to the 
“p-wave” spin density of conduction electrons, inducing a topologically non-trivial insulating ground- 
state. Based on the analysis of the charge- and spin-excitation gaps, the string order parameter, and 
the spin profile in the groundstate, we show that, at half-filling and low energies, the system is in the 
Haldane phase and hosts topologically protected spin-1/2 end-states. Beyond its intrinsic interest 
as a useful “toy-model” to understand the effects of strong correlations on topological insulators, 
we show that the p-wave Kondo-Heisenberg model can be implemented in p—band optical lattices 
loaded with ultra-cold Fermi gases. 


PACS numbers: 73.20.-r, 75.30.Mb, 73.20.Hb, 71.10.Pm 


Introduction. Topological Kondo insulators (TKI) are 
a type of recently proposed materials where strong in¬ 
teractions and topology naturally coexist m- Within 
a mean-field picture ISi], TKIs can be understood as a 
strongly renormalized /-electron band lying close to the 
Fermi level, and hybridizing with the conduction-electron 
d—bands. At half-filling, an insulating state appears due 
to the opening of a low-temperature hybridization gap 
induced by interactions at the Fermi energy. Due to the 
opposite parities of the states being hybridized, a topo¬ 
logically non-trivial ground state emerges, characterized 
by an insulating gap in the bulk and conducting Dirac 
states at the surface. At present, TKI materials, among 
which samarium hexaboride (SmBg) is the best known 
example, are under intense investigation both theoreti¬ 
cally and experimentally i7Hin] 

In order to gain further intuition into the effect of 
strong interactions, recently Alexandrov and Coleman 
in proposed an analytically tractable model for a 
one-dimensional (ID) TKI, i.e., the “p-wave” Kondo- 
Heisenberg model (p-KHM), consisting of a chain of spin- 
1/2 magnetic impurities interacting with a half-filled one¬ 
dimensional electron gas through a Kondo exchange [see 
Fig. Ha)]. The peculiarity of this model, which makes it 
crucially different from other one-dimensional Kondo lat¬ 
tice models studied previously |12fE5] , is that the Kondo 
exchange couples to the “p-wave” conduction-electron 
spin density, allowing for effective next-nearest neighbor 
hopping processes in the conduction band accompanied 
by a spin-flip [see Fig. [Il:a)[. Using a standard mean- 
field description Hi], the above authors found a topo¬ 
logically non-trivial insulating groundstate (i.e., a class- 
D insulator |24H2t)| ) which hosts magnetic states at the 
open ends of the chain. Soon after, two of us studied this 
system using the Abelian bosonization approach com¬ 
bined with a perturbative renormalization group analy¬ 



Figure 1: (a) Representation of the p—wave Kondo- 

Heisenberg model (p-KHM) in one dimension. The upper 
leg represents the conduction electron p-band, and the lower 
leg corresponds to a spin-1/2 Heisenberg chain. The Kondo 
exchange Jk couples a spin Sj with the p-wave spin density 
in the conduction band [see Eq. Q). (b) Microscopic model 
effectively realizing the p-KHM at low energies, and allowing 
an experimental implementation in p—band optical lattices. 
At half-filling, the fermionic sites in the lower-leg behave as 
spins due to a strong on-site Hubbard repulsion U. The direct 
hopping across a given rung vanishes due to the different par¬ 
ities of the orbitals. The unusual p—wave Kondo interaction 
in[^a) originates in second-order hopping processes in V, i.e., 
Jk = jU (see Appendix [a| . 

sis, revealing an unexpected connection to the Haldane 
phase at low temperatures m- The Haldane phase is a 
paradigmatic example of a strongly interacting topolog¬ 
ical system, with unique features such as topologically 
protected spin-1/2 end-states, non-vanishing string order 
parameter, and the breaking of a discrete Z 2 x Z 2 hid¬ 
den symmetry in the groundstate [HHSD]. The striking 
results in Ref. m indicate that ID TKI systems might 
be much more complex and richer than expected with the 
naive mean-field approach, and suggest that they must 
be reconsidered from the more general perspective of in¬ 
teracting symmetry-protected topological phases EBisa. 

In this Letter we study the groundstate properties of 
the finite-length p-KHM in one dimension using the Den- 
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sity Matrix Renormalization Group (DMRG) |33[ IM] , 
Our results indicate that the system is a Haldane insu¬ 
lator with protected spin-1/2 end-states and finite string 
order parameter, therefore supporting the predictions of 
Ref. [57]. We also propose that this exotic model could 
be realized in p—band optical lattices loaded with ultra¬ 
cold Fermi gases, which would allow for controlled exper¬ 
imental studies of TKIs in the lab. 

Model. The Hamiltonian of the p-KHM is H — 
Hi + H 2 + Hk [TTl I57j . where the conduction band 
is represented by a L-site tight-binding chain Hi = 

-t Ylf=i,cr (pj,crPi+i,o- + H.c.^ with p] the creation op¬ 
erator of an electron with spin a at site j with spatial 
p—symmetry [upper leg in Fig. [^a)]. The Hamiltonian 
H 2 = Jh J2jZi Sj ■ Sj-i-i [bottom leg in Fig. |^a)[ cor¬ 
responds to a spin 1/2 Heisenberg chain, and Hx is the 
Kondo exchange coupling between Hi and H2 mi HZ] 

L 

Hk = ( 1 ) 

i=i 


with Jk > 0. This Kondo interaction is unusual in that 
it couples the spin Sj in the Heisenberg chain to the 
“p-wave” spin density in the fermionic chain at site j, de¬ 


fined as TTj = J 2 a,l 3 { 


p]+i,, 


-pI- 


j-i,< 




'J - ^a,p\ ^ 

where po.o- = PL+i.a = 0 is implied, and where 
(To,p is the vector of Pauli matrices. Eq. Q can 
be written as Hk = H^k^ + where IPk'^ = 

^ Sj ■ (sj-i + Sj+i) contains the coupling of a lo¬ 
calized spin with the usual spin-density at site j ± 1 in the 
conduction band (where Sj = pPj 

H 


(2) _ _Jk 




'^a,/3 Pj + l,a{^^ )Pj-i,P + H.C. 
describes a different type of processes characterized by a 
non-local (i.e., next-nearest neighbor) hopping accompa¬ 
nied by a spin-flip. 

We have studied the groundstate properties of H by 
means of DMRG. In our implementation we have kept 
m = 800 states, which allowed to achieve truncation er¬ 
rors in the density matrix of the order of 10“^^. The 
DMRG method has been used previously to describe the 
standard ID Kondo lattice model at half-filling HMH], 
where a topologically trivial, fully gapped groundstate 
was obtained. For the p-KHM, where a topological in¬ 
sulator groundstate was predicted miiiz], there are no 
DMRG studies to the best of our knowledge. Intuitively, 
we expect that the charge and spin gaps in this model 
vanish in the limit Jk —>■ 0, as the Hamiltonians Hi and 
H 2 are separately gapless in the thermodynamic limit. 
According to the bosonization analysis in the limit of 
small Jk, both gaps are favored when the velocities of 
the gapless spinon excitations described by Hi and H 2 
coincide m- Intuitively, the term Hk becomes more ef¬ 
fective to couple the spin degrees of freedom in iJi and 
H 2 when they fluctuate coherently (i.e., same spinon ve¬ 


locities). The spinon velocity in the conduction band is 
equal to the tight-binding Fermi velocity vi = vp = 2t, 
and in the Heisenberg chain is V 2 = ttJh/2 |35I [36] . 
and therefore we conclude that the optimal situation 
in order to maximize the effect of Hk corresponds to 
Jh = 4t/7r « 1.27 t, which we choose in all our subse¬ 
quent calculations. In what follows, we characterize the 
groundstate by analyzing the charge and spin gaps, the 
string order parameter, and spin profile along the chain. 

Charge gap. Spin-flip scattering generated upon in¬ 
creasing Jk induces gapped charge- and spin-excitations 
in the system at half filling im HZ]. Although these 
gaps are not direct evidence of the topological nature 
of the groundstate, their study is important to charac¬ 
terize the p—KHM insulating phase and to test the pre¬ 
dictions of bosonization |5Z]. Using the hidden SU (2) 
charge pseudo-spin symmetry of the model at half-filling, 
we can compute the charge gap of a L—supersite sys¬ 
tem as Ac (L) = {N = L+ 2)- {N = L) 

dzids]. Here, a “supersite” j refers to the combination of 
a spin Sj and the fermionic site in each rung. is the 
z-projection of the total spin in the system, computed as 
where T^- = Sj+Sj. Finally, E^^ (N) 
is the groundstate energy of a system with N electrons 
in the conduction band, and projection M^. While pre¬ 
vious results on the standard ID Kondo lattice predicted 
a linear dependence Ac oc Jk HU [IS], here the presence 
of the Heisenberg coupling Jh changes this behavior as 
it cancels the first order contribution m- The leading 
second-order contribution to Ac can be physically un¬ 
derstood integrating out the “fast” spin fluctuations in 
the Heisenberg chain, which generate an effective repul¬ 
sive four-fermion interaction U' oc (Si.Sj )in the 
conduction p-band. This effective interaction produce 
umklapp processes which open a Mott insulating gap at 
half-filling [36]. In Fig. i we show the charge gap as a 
function of Jk ■ The system presents important finite-size 
effects in the limit of small Jk, and we therefore analyze 
our results with the scaling law Ac (L) « Ac (oo)-|-/3cT“^ 
in the case of large Jk {Jk/^ > 0.3) [HKH], whereas in 
the regime of smaller Jx/t < 0.3 the fits improve with 
the scaling law Ac (L) « -^A^ (cjo) -|- l3cL~'^ (see inset in 
Fig. §. This fitting procedure allows to extract Ac ( 00 ), 
the value of the charge gap in the thermodynamic limit, 
as a function of Jk (see Fig. |^. The solid (red) line 
is a quadratic law Ac (cxd) = acJx which fits the data 
reasonably well at small Jk, confirming the dependence 
predicted by bosonization m- 

Spin gap and spin-1/2 end states. We now focus on 
the spin degrees of freedom, where the p-KHM has the 
most interesting properties. Intuitively, the physics of 
the problem can be simply understood: the antiferromag¬ 
netic Kondo exchange along the diagonal rungs effectively 
forces the spins to align ferromagnetically across the 
rungs, even in the absence of a direct coupling |27| . This 
situation favors the formation of a local triplet in each 
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Figure 2: Charge gap Ac vs Jk- The circles are the val¬ 
ues of the gap in the thermodynamical limit Ac (oo), ob¬ 
tained after finite-size scaling (see inset). The solid (red) 
line is a fit Ac (oo) = UcJkj valid for small Jk, based on 
the bosonization analysis of Ref. m- The dashed lines are a 
guide to the eye. Inset: Finite-size scaling using the scaling 
laws Ac (L) « Ac (oo) -I- /3cT~^ for Jk/I > 0.3 [161 [19], and 
Ac (L) « (oo) + /3^L-2 for Jk/I < 0.3. 


supersite, and the system mimics the properties of the 
spin-1 Heisenberg chain m or the ferromagnetic Kondo 
lattice model |15l [38] , which are examples of systems re¬ 
alizing an insulating Haldane groundstate. A hallmark of 
this phase is the presence of two topologically protected 
spin-1/2 magnetic states at the ends of the chain (i.e., 

lt)L ® lt)_R) lt)L ® ® lt)_Ri and 1/)^ (g) 1/)^ ), 

which arrange into degenerate triplet and singlet linear 
combinations. As a result, the first spin-excitation gap 
(L) = {N = L)- {N = L), tends to 

zero for L —> oo. For a finite-L chain, however, the over¬ 
lap of the end-states wavefunctions removes this degener¬ 
acy exponentially as (L) oc where ^ oc 

is the localization length for the magnetic end-states, 
and the groundstate for N even (odd) corresponds to 
the singlet S' = 0 (triplet S = 1) [34] . In our case, at 
small Jk the localization length ^ becomes of the or¬ 
der of the system size (^ ~ L), and it was not possible 
to obtain a conclusive scaling behavior for even 

for the largest systems we have simulated (L = 120). 
On the other hand, the gap (same definition as 

above changing = 1 —)■ = 2) can be identified 

with the Haldane gap of the system, and physically in¬ 
volves spin excitations which live in the bulk (see Fig. 
1^. In this case, the scaling analysis is simpler as it is 
free from edge effects, and we have used the scaling law 
A, (L) « (oo) -I- /3sF“2 for all values of Jk (see bot¬ 
tom inset in Fig. |^. The values of Ag (oo) are shown 
in Fig. In contrast to the case of the charge gap, 
here the analytic dependence of Ag (oo) on the parameter 
Jk is technically more challenging to obtain within the 


Figure 3: Spin gap Aa^’°^ (physically corresponding to the 
Haldane gap in the system), as a function of Jk- The ex¬ 
trapolated spin gap Aa^’°^ (oo) is obtained after a finite-size 
scaling analysis (see bottom inset). The dashed lines are 
a guide to the eye. Top inset: Spatial profile of (TJ) = 
i.e., the z component of the spin in the 
supersite j, computed with the groundstate of the subspace 
with total = 1 for L = 80. The presence of the topo¬ 
logically protected spin-1/2 states at the ends of the chain 
is clearly seen. Botto m inset: Finite-s ize scaling using the 
scaling law Aa (L) Ri -^/Af (oo) + PsL~^. 


bosonization formalism, and is beyond the scope of this 
work. Nevertheless, our numerical results yield a power- 
law dependence Ag (oo) oc J)/, with exponent v = 2. 


We next investigate the presence of topologically pro¬ 
tected spin-1/2 end-states which, as mentioned before, 
is a crucial feature of the open Haldane chain. In the 
upper inset of Fig. we show a spatial profile of the 


0 —projection of T^, i.e., (T?) = | T, 




M^ = l 

g 


where is the groundstate with total spin = 

1, for Jic/t = 1 and Jiy/t = 2. For these large values 
of Jk (which are beyond the validity of the bosonization 
analysis 1221) the end-states are clearly visible and show 
a small localization length a fact that prevents them 
from overlapping, producing negligible finite-size effects. 
Since we are working in the subspace = 1, and since 
the spin profile is symmetric under space inversion, we 
conclude that the spin accumulation at each end is 1/2, 
corresponding to the configuration where the topological 
spin-1/2 end-states is It)/, ® It)7? 1^ . 


String order parameter. The most fundamental sig¬ 
nature of the Haldane phase is, however, the emer¬ 
gence of a finite string order parameter |39| . a quantity 
deeply connected to a broken hidden Z 2 x Z 2 symme¬ 
try m- This quantity is a smoking-gun for the pres¬ 
ence of the Haldane phase, and therefore is the most im¬ 
portant for our present purposes. Using the above def¬ 
inition of Tj, the string order parameter is defined as 
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Figure 4: String order parameter vs Jk- Throughout 

the whole studied regime, remains hnite, indicating 

the presence of a Haldane-insulating phase. Inset: Spatial 
dependence of Ogtring (d-) vs the distance d — |1 — m|, where 
I and m have been taken symmetrically about the center of 
the chain. 


C’s“tring (^ - H Due to the 

SU (2) spin-symmetry of the model, it is enough to cal¬ 
culate the computationally simpler component a = z. 
We have computed Oftring ~ taking the sites I and 
m symmetrically about the center of the system in order 
to minimize the effect of the edges. Note in the inset of 
Fig. I^that 0|tring converges rapidly as a function of 
the distance d =\l — m\io the ID-bulk value In 

the main Fig. 4 we show vs Jk, which remains fi¬ 

nite throughout the whole studied regime. This indicates 
the presence of a Haldane-insulating phase even beyond 
the regime of small Jk where the bosonization analy¬ 
sis in Ref. m is valid. This result, together with the 
confirmation of the presence of spin-1/2 end states, are 
the most important results of this paper, as they provide 
conclusive evidence that the p-KHM realizes a Haldane 
insulating phase. 

Realization in experimental systems. Recent ex¬ 
perimental [4fil - l42] and theoretical |43l |44] works on 
optical lattices with higher orbital p—bands suggest 
that the p-KHM could be realized in the laboratory. 
More specifically, the tight-binding Hamiltonian Hi 
could be simulated populating the first excited en¬ 
ergy level with p—symmetry at each site in a ID 
optical lattice (see upper leg in Fig. Bb)). The 
Heisenberg chain H 2 could be simulated with a half- 
filled Hubbard model in the Mott insulating phase 

US SB], i.e., i?Hubbard = + H.c.^ -b 

~ 5 ) ~ 5 ) ™ ^be limit U :$> t' (see 

lower leg in Fig. [Ub)), where U is the on-site Coulomb 
interaction which could be tuned via a Feschbach reso¬ 


nance. Here ^ is a creation operator of a fermion with 

spin a at the s—orbital site j, and n/ / = s/ „Sj^cr the 
fermion occupation. The unusual Kondo exchange Hk 
would naturally arise in this situation if a microscopic 
single-particle hopping between Hi and i?Hubbard is al¬ 
lowed [see Fig. Bb)]. Due to the different parities of the 
orbitals involved, the matrix element connecting the sites 
on the same rung (i.e., along the vertical direction) van¬ 
ishes. The leading contribution therefore corresponds to 
the matrix element V coupling s and p orbitals along a di¬ 
agonal rung, i.e., {Vj+i,<y - Pj-i,a) + 

H.C., where the crucial sign inside the parentheses is a di¬ 
rect consequence of the p—wave nature of the conduction 
band states. The equivalence between the more “phys¬ 
ical” Hamiltonian H' = Hjiubbard + ^^s-p + Hi and the 
p-KHM can be rigorously shown by the means of a canon¬ 
ical (i.e., a generalized Schrieffer-Wolff) transformation 
T = e**^, where the operator S = S {t' ,V) must be cho¬ 
sen so as to eliminate first order contributions in t' and 
V. The procedure is standard and here we only out¬ 
line the main steps (see Appendix [A| for details). As¬ 
suming the limit t',V <C U, we can expand the expo¬ 
nential in H and truncate the series at second order in 
(t'/U) and (V/U), therefore obtaining H = H'T « 
H' H'] -b ^ [iS, [iS, H']]. We choose the transforma¬ 
tion to be 5 = -i [(iJ/- - iJjT) -b 2 {H+p - iy^lp)] /U', 
with H+ = -t' 


N (s) t 




1 - n. 


(D 


(with 


(Pi-l-l.cr Pi —I,a) ^1 

the notation (ij) indicating that i and j are nearest- 
neighbor sites), and where H/i = {H^) and H~p = 
(Ffgtp)^- It is easy to check that the first order contribu¬ 
tions cancel, and therefore H 2 = {Hi[Hi^^ V and 

Hk = —yR (H^pH^p) V, where V is the projector onto 
the lowest subspace of iFnubbard (see Appendix . The 
connection between both models is completed identifying 
the parameters as Jh = 4F^/[/ and Jk = 8R^/C/. 


Note that this proposal is different from other theo¬ 
retical proposals to simulate the standard Kondo lattice 
model in ID optical lattices |47l |48| . 

Conclusions. We have studied the p-KHM, a theoreti¬ 
cal “toy-model” introduced to describe a one-dimensional 
topological Kondo insulator, by the means of DMRG, 
and have calculated various quantities characterizing the 
properties of the groundstate at half-filling. We have 
shown strong numerical evidence (based on the analysis 
of the charge and spin gaps, the spin profile, and the 
string order parameter) that the p-KHM realizes a Hal¬ 
dane insulating phase at low temperatures, as predicted 
in Ref. [27]. Our results indicate that the topological 
properties of this model fall beyond the scope of the 
non-interacting topological classification |24bl26| . which 
is unable to reveal the true topological structure of the 
groundstate. Finally, we have proposed that the unusual 
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p—wave nature of the Kondo interaction could be physi¬ 
cally realized in experiments with ultracold Fermi gases 
loaded in p—band optical lattices. 
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Appendix A: Derivation of the p-KHM by a Canonical Transformation 


In this Appendix we provide a derivation of the p-KHM Hamiltonian H in the main text by the means of a canonical 
transformation. To that end, we start from the microscopic Hamiltonian H', consisting of a fermionic Hubbard ladder 
with s and p orbitals along the legs, and depicted in Figure 1(b) in the main text: 


■^Hubbard “1" ^s-p 

(Al) 

L-1 


— t + + H.c.^ 

(A2) 


(A3) 

L-1 


{p].aP3+i,<y + H.c.) 

(A4) 

J —1,(T 


(P3+h<T - Pj-i,<y) + H-c., 

(A5) 


J.o- 


Note that the system has electron-hole symmetry. Here, s] ^ creates a fermion with spin projection cr at site j in the 

Hubbard leg and = sj ^Sj^a is the corresponding fermion-number operator. The operator ^ creates a fermion 
with spin cr at site j in the p—orbital conduction band, represented by a simple tight-binding model Hi. The term 
Hs-p couples the two fermionic legs, and due to the symmetry properties of the s— and p—orbitals, the direct hopping 
across the rungs is zero. Therefore, the most important hopping process occurs between a fermion Sj^a- and the linear 
superposition with p—wave symmetry cx (pj+i^g- — Pj-i^^) in the conduction band. 

The idea is to derive an effective low-energy model in the limit U {t', V}. To that end, we split the Hamiltonian 
H' into 


H' = Hf + H,.p + Hu + Hi, (A6) 

where 

L 

Hr = -t' (sl^Sj,a + H.c.^ , 

Hu = uY - 0 

The first two terms in ( |A6| ) will be considered as perturbations to Hu, in the regime {f, V} U. 

We now start from the atomic limit in the Hubbard leg, i.e., t' = V = 0, and identify the atomic singly-occupied 
states \aj) = |0) =t)i) as forming the lowest-energy subspace at site j, while the |0j) (empty) and \dj) = 

"ps] ^ |0) (doubly-occupied) form the excited subspace. We now introduce projectors onto each of the 4 atomic states: 


(A7) 

(A8) 


n.»=(i-"S)(i-";3). (“) 

= "S3"i3. (A“) 

Pj.r = »;•> (l - n'5) , (All) 

Pu = n<:> (l - nS:') . (A12) 

Note that while all projectors commute with Hu, the kinetic terms Hr and Ffg-p cause transitions among subspaces. 
Using that Ij = can write the kinetic terms as Hr = -I- H^ + H^,, 

and Hs.p = i7s-p (Ej ./3 = H+p + iJ^p + where 
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N 


= -t' (l - + <^s\ „s. 


(s) t 


N 


(i - 


^t' =-*' Y [(l “ 4,<rS^,<ynl"l + (l - slaSj,an^"l 

{ij),<^ 

= (p^+l,<y - Pi-1, a) + (pl+l,a - PLl,cr) Si,a (l - 

2,(7 

H^-P = ^Y [(^i+l.o- “ pI-1.o) Si,an^il + (l “ ”-i!ff) ®i,cr (Pi+l.o- “ Pi-l,<j) 


(A13) 

(A14) 

(A15) 

(A16) 


Physically, the term with supraindex “+” produce transitions from the lowest subspace to the excited subspace, while 
those with “—” restore excited states to the lowest subspace. On the other hand, the terms labelled with “0” do not 
change the subspace, and since we assume a half-filled conduction band, they will identically vanish and it is not 
necessary to write them explicitly here. We now note the following important relations 


H^ = {H+)\ (A17) 

H-^={H+^)\ (A18) 


which will be useful in what follows. 

We now introduce a canonical transformation in Eq. 
taneously get rid of the terms at first order in t' and V: 


(All, such that in the transformed representation we simul- 


H = 

= H' + t[S,H'] + '^[S,[S,H']] + ... 


(A19) 

(A20) 


We want to choose S in such a way that H does not connect different Hubbard subbands. Note that this cannot be 
achieved at infinite order in the expansion in powers of S in Eq. (A20l, but we will be content if we can eliminate 
the contributions at order O {t') and O {V) that mix the subbands. We now write the expansion in Eq. (A20l in the 
more suggestive form 


H = H+ + + i7-p + * [5, Hu] 

+ Hi+ Hu + i [5, H^ + H^,] +i [5, H^_^ + iJg_p] 

+ — [5, [5, Hu\\ + (other less important terms) 


(A21) 

(A22) 

(A23) 


We will require that the first line (A21l in the above equation vanishes. It is then clear that S must be O {f /U) 
O {VIU). Using the following results 


(s) t 





{p\+Ua - P\-1,^ (l - , 




(l - ng) , Up 


2 / \ 


1 

--TL'- - 


1 

2 ’ 


-p\P\,a (Pi+1.<T -Pi-l,cj) , 


^ {p\+Ua - P\-I,a) S^,a (l - UPj , 


(A24) 

(A25) 

(A26) 

(A27) 


it is easy to check that 
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(A28) 
(A29) 

Then, it follows that the choice 


[H^,Hu] =tUH^, 


2i 


S = -jy{H+-H^)--{H+p-H- 


(A30) 


exactly cancels line (A21l. 

The relevant part of the Hamiltonian at low energies is then obtained projecting H onto the lowest Hubbard 
subband. This is formally done applying the projector V = (^i,t + )> which eliminates certain terms in Eqs. 

(A22l and (A23). The resulting effective Hamiltonian at lowest order in t/U and V/U is therefore 


H = H,+v\Hu + i[S,H+ + H^] +1 [S, H+p + i7-p] + ^ [5, [5, Hu]] \ V 


= H,+Hu-^V {H-H+) V-^V (H-pH+p) V. 


We now replace the expressions for and H^ [Eqs. (A13l-(A16)] into the above equation and obtain 


V {^H^, V — {t ) ^ ^ \Vi,a'Pi+l^a + T’i,a'Pi+l,(T 

1 , (7 

V {H-pH+p) V = 2V^Y. ^ E 

2 , (7 2,(7 

Pi_i {Pi+l,a Pi—l,a) T (Pi+l,a Pi—l,a^ iPi+l,tT Pi—1 ,ct) 

Using that + Vi^^Pi+i^^ = —2SfSfj^i + n^®^n-^\/2, and the Schwinger-fermion representation 


(A31) 


(A32) 


(A33) 


= 


('*) (s) 


5+ = 


2 

t 


is a faithful representation of a spin-1/2 operator, we can write the effective Hamiltonian as 


H = Hu + H, + ^-^f;^ 


0.0. , st+,s-+ s-^,st 1 

Oi -h 


8 U2 ^ ^ - p\_^^^) (Pz-Hi.t - P*-i.t) - (pI+1,; - pLi. J - P*-id) 


8U2 

+ irE 


Sf (pl+i,i -pLi.i) (P*-si.T -K-i.t) ^ S, (pI+ 1^^ -pLi,t) (Pi-si,l -Pi-i.l) 


where we have neglected the constant ^ Defining the effective parameters 




Jk = 


U ’ 
8U^ 
“cT’ 


we note that this Hamiltonian corresponds to the p-KHM considered by Alexandrov and Coleman m- 


(A34) 

(A35) 

(A36) 


(A37) 


(A38) 

(A39) 





















